arXiv: 1508.0555Ivl [astro-ph.HE] 23 Aug 2015 


PROCEEDINGS 

SCIENCE 



VERITAS detection of y-ray flaring activity from the 
BL Lac object 1ES 1727+502 during bright 
moonlight observations 


Matteo Cerruti* for the VERITAS Collaboration 

Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA 
E-mail: [natteo . cerruti@cfa . harvard. ed~ 

During May 2013, a gamma-ray flare from the BL Lac object lES 1727-\-502 (z=0.055) has been 
detected with the VERITAS Cherenkov telescopes. This detection represents the first evidence 
of very-high-energy (E>100 GeV) variability from this blazar and has been achieved using a 
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gral flux is about five times higher than the archival VHE flux measured by MAGIC. The detec¬ 
tion triggered additional VERITAS observations during standard dark-time and multiwavelength 
observations from infrared to X-rays with the FLWO 48" telescope and the Swift satellite. The 
results from this campaign are presented and used to produce the first spectral energy distribution 
of this object during gamma-ray flaring activity. The spectral energy distribution is then fit with 
a standard synchrotron-self-Compton model, placing constraints on the properties of the emitting 
region in the blazar. 
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1. Introduction 

Active galactic nuclei (AGN) of the blazar class represent the most common source of very- 
high-energy (VHE; E > 100 GeV) y-ray photons in the extra-galactic sky. Their emission is charac¬ 
terized by a broad spectral energy distribution (SED) from radio to y-rays, extreme variability and a 
high degree of polarization. In the framework of the unified AGN model they are considered AGN 
whose relativistic jet is aligned with the line of sight. The study of VHE blazars is complicated by 
their broadband emission and their rapid variability. Strictly simultaneous multiwavelength (MWE) 
campaigns are the key to study blazar physics. 

The cameras of imaging atmospheric Cherenkov telescopes are made of photomultiplier tubes 
(PMTs), which cannot safely operate in bright moonlight conditions. This constraint particularly 
affects the studies of blazars (and other variable sources), limiting the organization of MWE cam¬ 
paigns, or prohibiting follow-up observations of a flare. Among the current generation of lACTs, 
only MAGIC and VERITAS observe under moderate moonlight (<35% moon illumina¬ 

tion). FACT, whose camera is instead composed of solid-state Geiger-mode avalanche photodiodes 
(also called silicon-photomultipliers) 0, is also capable of performing moonlight observations. 

In 2012, the VERITAS collaboration began a new observing program under bright moonlight 
by reducing the high voltage (RHV) applied to the PMTs, or by utilizing UV-transparent filters. 
The details of the observing strategy are discussed in [p^]. In this proceeding we present the detec¬ 
tion of VHE flaring activity from the blazar lES 1727+502 during May 2013, at a flux of roughly 
five fimes fhe archival VHE flux measured by MAGIC [^. The high-flux sfafe was inifially defected 
during brighf moonlighf observations, which represenfs an innovafion for VERITAS. 

The blazar lES 1727+502 is a nearby high-frequency-peaked BE Eac objecf (z = 0.055), dis¬ 
covered as a y-ray source by Fermi and as a VHE source by MAGIC [^]. The MAGIC collaboration 
reported an integral flux of (5.5 ± 1.4) x 10^^^ cm^^ s^^ above 150 GeV and a specfral index of 
2.7±0.5 ||]. 

Eurfher defails on fhe VERITAS mulfi-wavelengfh campaign are presenfed in [^. 


2. VERITAS instrument 

VERITAS is an array of four lACTs located af fhe Whipple Observatory in soufhern Arizona 
af an alfifude of 1.3 km above sea level. Each telescope is of Davies-Coffon design wifh a 12-m 
diameter reflecfor, and fhe array is arranged in a diamond wifh ~ 100 m fo a side. Each VERITAS 
felescope is insfrumenfed wifh a camera made up of 499 PMTs each, wifh a lofal field-of-view of 
3.5°. VERITAS is able fo reconsfrucf y-ray emission from abouf 85 GeV fo >30 TeV, and a single- 
evenf angular resolution of 0.1° af 1 TeV. During sfandard observations (i.e. when fhe moon is 
< 35% illuminafed) a source wifh an infegrafed flux of 1% of fhe Crab Nebula flux can be defecfed 
af fhe 5a level in ~ 25 hours, and a 5% Crab source in less fhan 2 hours. More informafion on fhe 
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VERITAS array can be found in Q, p^ . 

In the RHV observation mode, the PMT voltages are reduced to 81% of their standard values 
during dark-sky observations, allowing VERITAS to operate when the moon is 35-65% illuminated. 
The sensitivity of the RHV observation mode is similar to the standard VERITAS sensitivity, albeit 
with a higher energy threshold (~ 200 GeV). The cuts used in this RHV analysis would allow for 
a 5% Crab source to be detected in less than 2 hours. 


3. VERITAS observations of lES 1727+502 

The VERITAS observations that allowed for the detection of lES 1727+502 were taken be¬ 
tween May 1, 2013 and May 7, 2013. There were additional data taken on May 18, 2013 that did 
not result in a detection. After quality selection, approximately 6 hours of data remain. Of these, 
3 hours were taken in RHV mode on the first two nights of the exposure. All observations were 
made in “wobble” mode, wherein the telescopes were pointed 0.5° away from the target to allow 
simultaneous measurements of the target and background regions. All data in the 2013 dataset 
were taken at elevations between 65° and 70°. 


The analysis of the complete 2013 dataset used the “reflected-region” background model 0 
resulting in 159 ON events and 850 OFF events with a background normalization factor a of 0.077, 
yielding a detection significance of 9.3a []A|. The reconstructed VHE spectrum is shown in Fig. |^. 
In the same figure we show data points corrected for absorption by extragalactic background light 
(EBL) using the model in [^] (for z=0.055). Data are well fitted with a power-law function be¬ 
tween 0.25 TeV and 1.6 TeV; the /NDE is 1.66/2 (P 44%). The observed spectrum is given 




= (7X±1.1)x10-'Mto) 


-2.1±0.3 


cm 


-Ir- 


TeV ^ All errors in the aforementioned fits 


620 GeV y 

are statistical only. The cumulative systematic errors on the flux normalization and spectral index 
are conservatively estimated to be 30% and ±0.3, respectively. 


The light curve for the VERITAS observations is shown in the top panel of Eig. The peak in¬ 
tegrated flux above 250 GeV is (1.6 ±0.4) x 10~'^cm^^s^\ which corresponds to 9.5% of the Crab 
Nebula flux. The upper limit is at the 95% confidence level and represents 9.4 x 10^'^cm^^s^^ 
(5.6% Crab). The light curve between May 01 and May 07 (MJD 56413-56419) can be fitted with 
a constant resulting in a flux of (1.1 ±0.2) x 10^^'cm^^s^^ (6.3% Crab) and a x}/AIDE = 4.73/3 
(P ~ 19%). 


The integral flux above 250 GeV corresponds to about five times the flux measured by MAGIC, 
and represents the first evidence of VHE variability in lES 1727+502. The VERITAS observations 
between May 01 and May 07 are consistent with a constant flux. However, when including the 
measurement on May 18 in the fit, the resulting x'^/NDE is 12.0/4 (P ~ 1.8%). Thus, a constant 
flux is excluded at the 2.4a level, indicating that the flare may have ended at some point after the 
last VERITAS detection on May 07. 
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Figure 1: Spectrum for the complete lES 1727+502 dataset. Power-law fits to both the uncorrected (black) 
and EBL-corrected (red) points are provided. The blue shaded region represents the result from MAG/C [3]. 
For further details see Q| . 


4. Multiwavelength observations of lES 1727+502 

Following the detection of 7 -ray flaring activity from lES 1727+502, we triggered observa¬ 
tions at lower wavelengths, to characterize the SED of the source during the active state. The Swift 
satellite observed lES 1727+502 on May 7, 2013 and May 8 , 2013, providing spectral reconstruc¬ 
tion in soft X-rays (using the XRT telescope) and optical/UV (using the UVOT telescope). The 
ELWO 48" optical telescope observed the source on May 18, 2013, ten days after the end of the 
VERITAS flare. The MeV-GeV part of the SED is covered by the Eenni-LXT telescope, thanks to 
its monitoring capabilities. 


We performed two different analyses of Een72/-EAT data: we first studied the long-term av¬ 
erage spectrum of lES 1727+502, and then focused on the MeV-GeV emission spanning the in¬ 
terval of 7 -ray activity detected with the VERITAS telescopes. Using data from August 4, 2008 
to August 1, 2013, lES 1727+502 is detected with a significance of about 14 standard deviations 
over the background, and its emission between 100 MeV and 300 GeV can be parametrized by 
a power-law function with index Eoj-sooGeV = 1-91 ±0.08 and differential flux ‘Eo.i-sooceV = 
(2.1 ±0.2) X 10^^^ cm^^ s^^ MeV^', estimated at the decorrelation energy £'o;o.i- 3 ooGey = 2136 
MeV. This value is consistent with the average spectrum provided in the 2EGE catalog [IE]. A sec¬ 
ond analysis was performed in order to produce a measurement in the GeV energy band spanning 
the same interval as VERITAS observations. Only Eenni-LKT observations taken between May 1 
and May 7, 2013 are considered. lES 1727+502 is not detected by Eenni-EAT in this very short 
period and only upper limits on its GeV emission can be computed. 


The detection of the high-flux state at VHE by VERITAS triggered X-ray and UV observa¬ 
tions by the Swift satellite, which observed lES 1727+502 on May 7 and May 8 , 2013 for a total 
live time of 6.8 ks. The XRT light curve above 0.3 keV (in counts per second), corrected for the 
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Figure 2: Light curve of lES 1727+502 during the May 2013 multiwavelength campaign. From top to 
bottom; the VERITAS light-curve, the 5w!/f-XRT light-curve and the 5w!/f-UVOT / FLWO 48" light-curve. 
See for details. 


exposure and for background, is shown in Fig. ^ No significant variability is detected within a 
single observation, nor between the two observations. The average count rate measured by XRT is 
0.83 ±0.12 counts per second. Given the lack of variability, and in order to improve the statistics, 
the data are summed and rebinned assuming a minimum of 50 counts per bin. Data below 0.3 keV 
are excluded, and the last significant spectral bin extends up to 8 keV. The average XRT spectrum 
is well fitted by an absorbed (Nh = 2.75 x 10^° cm^^, as measured by 0) broken-power-law func¬ 
tion. The best-fit parameter values are Fi = 2.01 ^q }q, r 2 = 2.44^q}q, Etreak = 1-21^0 20 keV, and 
normalization C = ( 5 . 3 +Q 2 ) x 10^^ cm^^ s^^ keV^^ The x^/EfDF value is 86/84. 

The UVOT telescope, on board the Swift satellite, observed lES 1727+502 at optical and ul¬ 
traviolet wavelengths, simultaneously with XRT. All measurements were performed using the six 
available UVOT filters. As part of a long-term optical program of monitoring of VHE blazars, 
lES 1727+502 is regularly observed by the automatic ELWO 48" telescope, located near the VER¬ 
ITAS site. There are no observations performed simultaneously with VERITAS: the observation 
closest to the VERITAS detection of lES 1727+502 was performed on May 18, 2013, eleven days 
after the last detection by VERITAS using B, r’, and i’ filters. The Galactic extinction is taken into 
account assuming Eb-v = 0.037, a value consistent with the Nh value used for the X-ray analysis. 
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For the purpose of subtracting the host-galaxy contamination, we made use of the recent results 
from lllsl]. 


5. SED modeling 

The SED of lES 1727+502 is shown in Fig. In the context of the synchr'otron-self-Compton 
(SSC) model, the two components of the blazar SED are associated respectively with synchrotron 
emission from leptons (e^) and inverse-Compton scattering of the particles off synchrotron radia¬ 
tion from the same lepton population. The emitting region is a spherical blob of plasma (charac¬ 
terized by its radius R) in the relativistic jet, moving towards the observer with Doppler factor 5, 
and filled with a tangled, homogeneous magnetic field B. The parficle population is paramefrized 
by a broken power-law function, and if carries six free paramefers: fhe minimum, maximum, and 
break Eorenfz factors of fhe parficles fmax^ Ybreak)< the fwo indices (ai and a 2 ), and fhe nor- 
malizafion facfor K. The minimum Eorenfz facfor of fhe leptons can be fixed af a reasonably low 
value wifhouf affecfing fhe modeling. The remaining eighf free parameters can be constrained by 
observations, as discussed, for example, in 0 (and references therein). 

Given the lack of a simultaneous Femii-LPH detection, it is impossible to constrain the posi¬ 
tion and the luminosity of the inverse-Compton peak during the high state, and a unique solution 
for the SSC model cannot be provided. However, the synchrotron component is very well sampled, 
and can provide some constraints on the energy distribution of paiticles in the emitting region. At 
low energies, the subtraction of the host-galaxy contamination reveals the AGN nonthermal con¬ 
tinuum, which can be described by a power law from infrared to UV, with no sign of break. A fit 
of the FLWO 48” and UVOT data results in an index n\ = 1.64 ±0.09. In the case of synchrotron 
radiation, this index reflects directly the index of the underlying e^ population a\=2n\ — \, which 
is thus equal to 2.28 ±0.18. Similarly, the index of the X-ray power law below the X-ray break 
corresponds to a 2 = 3.0 ±0.2. 

The value of does not affect the modeling and it is therefore held fixed. However, values 
lower fhan 10^ would overesfimafe archival radio measuremenfs (see Figure]^, even fhough, given 
fheir nonsimulfaneify, fhese dafa should nof be considered as a sfrong consfrainf. An additional 
consfrainf on 7 ,„,„ is provided by fhe Fenni-LPVF nondefecfion: in fhe following we sfudy fwo dif- 
ferenf cases, for 7 ^, 7 , = 10^ and 5 x 10^. The value of jmax is consfrained by fhe break observed 
in fhe X-ray specfrum of IFS 1727+502. In fhe following, we express if as a function of ybreak'- 

Ymax 22 Ybreak- 


The number of free paramefers in fhe SSC model is five (5, B, R, Ybreak, and K), and wifh 
only four observables (fhe frequency and flux of fhe synchrofron peak, fhe VHE specfral index and 
fhe VHE flux af 620 GeV: v^yn-peak', vFy-syn-peak', ^VERITAS and vFv-yERiTAs), a unique solution 
cannof be provided. However, assuming a reasonable value of fhe Doppler facfor, if is possible fo 
sfudy fhe paramefer space of fhe remaining free paramefers. The besf-fif solution is compufed using 
fhe numerical algorifhm described in [^]. The Doppler facfor has been fixed fo 30. All fhe SSC 
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Figure 3: Quasi-simultaneous SED of lES 1727+502 (in blue; from low to high energies: FLWO 48", Swift- 
UVOT, Swift-XRT, Fermi-LAIT and VERITAS). Archival data (from the NED, Fermi-LAIT, and MAGIC) are 
included, and plotted in gray. The light-blue and pink curves represent the two sets of SSC models computed 
assuming 7 „,„, = 10^ and 5 x 10^. See |Q| for further details. 


models that correctly describe the lES 1727+502 SED are then recomputed, and plotted in Fig. 
^ The magnetic field is between 0.3 and 0.6 mG, with the emitting region size between 4x 10^^ 
and 7x10^^ cm, andjhr between 2x10^ and 3x10^. The emitting region is significantly out of 
equipartition {u^juB = lO^^'*), and the luminosity is of the order of lO'*"^ erg s^^ 


6. Conclusions 

We presented the results of a multiwavelength campaign carried out on the BL Lac object 
lES 1727+502 during May 2013, triggered by a VHE high-flux state detected by VERITAS. This 
represents the first detection of a blazar flare with VERITAS during bright moonlight observations. 

Within this campaign, no significant variability is detected at VHE, nor at lower energies (X- 
rays and optical). The VERITAS light curve is consistent with a constant flux (6.3% Crab above 
250 GeV, roughly five times the archival MAGIC detection); additional observations ten days af¬ 
ter the VERITAS campaign indicate (at a 2.4a level) that the high-flux state may have ended at 
some point after the last VERITAS detection. The quasi-simultaneous SED is fitted by a standard 
one-zone SSC model. Even though the nondetection by Eertni-EKT did not enable a full study of 
the parameter space, the measurements are fully consistent with particle acceleration by relativistic 
diffusive shocks and simple synchrotron and inverse Compton cooling, resulting in a power law 
with injection index around 2.2 and a spectral break of 1.0. 

In recent years the number of VHE blazars has significantly increased thanks to the current 
generation of Cherenkov telescope arrays such as VERITAS, MAGIC and H.E.S.S.. The broad¬ 
band emission and the rapid variability of blazars require prompt simultaneous multiwavelength 
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campaigns to fully understand the physics of their emitting region. The new VERITAS observing 
strategy under bright moonlight will be particularly useful for blazar science, significantly increas¬ 
ing blazar monitoring capabilities at VHE. 
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